A Microscopic View on the Mott transition in Chromium-doped V2O3 
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PACS numbers: 71.30.+h, 78.30.-j, 62.50.+p 

V2O3 is the prototype system for the Mott 
transition, one of the most fundamental phenom- 
ena of electronic correlation. Temperature, dop- 
ing or pressure induce a metal to insulator transi- 
tion (MIT) between a paramagnetic metal (PM) 
and a paramagnetic insulator (PI). This or re- 
lated MITs have a high technological potential, 
among others for intelligent windows and field 
effect transistors. However the spatial scale on 
which such transitions develop is not known in 
spite of their importance for research and appli- 
cations. Here we unveil for the first time the MIT 
in Cr-doped V2O3 with submicron lateral resolu- 
tion: with decreasing temperature, microscopic 
domains become metallic and coexist with an in- 
sulating background. This explains why the as- 
sociated PM phase is actually a poor metal. The 
phase separation can be associated with a ther- 
modynamic instability near the transition. This 
instability is reduced by pressure which drives a 
genuine Mott transition to an eventually homo- 
geneous metallic state. 

Hitherto, the phase diagram of V2O3, see FigQJa), was 
mainly established by means of transport measurements, 
i.e. detectingthe electronic response on a macroscopic 
spatial scale According to this phase diagram, V2O3 
exhibits a first-order MIT from a low-T antiferromag- 
netic insulating state (AFI) to a paramagnetic metallic 
state (PM) when T exceeds T N ~ 160 K. A different, 
paramagnetic insulating state (PI) is reached at T > Tjv 
if V is partially substituted by Cr. Within a narrow 
range of Cr content (0.005< x <0.018) a PI to PM tran- 
sition is obtained upon reducing the temperature across 
a Tmit > Tn with a slight discontinuity of the lattice pa- 
rameters within the same corundum phase. In particular 
(Vo.gsgCro. 011)20^ has attracted great interest since its 
discovery in 1969 [l|, because it provides the opportunity 
to span all phases, through two MITs. One (AFI o PM) 
shows both a change of crystal symmetry and of mag- 



netic state, the other (PM o PI) is iso-structural and a 
prototype of the correlation-driven Mott transition. 

It was soon evident that this simple scheme hides a 
more complicated scenario. In their conducting phase 
weakly Cr-doped samples show a poor-metal behavior 
(resistivity ~ 1(T 2 l]cm) 0,0]. Moreover the structure 
of (Vo.gsgCro. 011)203 is described as coexistence of an a- 
phase, with lattice parameters close to those of metal- 
lic V2O3 at 400 K, and a /3-phase with basically the 
same structure as at higher Cr content [3, 0]- hi ad- 
dition, EXAFS measurements [H suggest the essential 
role played by the local lattice strain which occurs in the 
Cr-doped compounds: Cr contracts the Cr-V bond, in- 
ducing a concomitant elongation of the V-V distances Q . 
As shown by theoretical calculation based on the Local- 
Density Approximation combined with Dynamical Mean 
Field Theory (LDA+DMFT) @], a longer V-V bond is 
associated with the PI insulating phase [3| . Besides these 
open issues concerning the T-induced MIT, further ques- 
tions are raised by studies under external pressure: in 
(Vo.gsgCro. 011)203, for instance, a conducting phase is 
obtained at room temperature if P > 3 kbar [7fl , and early 
investigations stated that applying pressure is equivalent 
to decreasing the Cr content [8| . Recently it was observed 
that different electronic structures for the metallic phase 
can be obtained when the PI-PM transition is induced 
by doping temperature or pressure ||, challenging the 
simple scenario of the McWhan phase diagram [lj. 

In this paper we explore the phase diagram of Cr- 
doped V2O3 from macroscopic to microscopic scales, by 
combining three spectroscopies with different lateral res- 
olution, i.e., Infrared (IR), Scanning Photoemission Mi- 
croscopy (SPEM), and X-ray Diffraction (XRD), with 
LDA+DMFT calculations. We find unambiguous evi- 
dence of an electronic phase separation in the T-induced 
PM phase. The poor-metallic behaviour is therefore ex- 
plained in terms of a coexistence of metallic and insu- 
lating domains stabilized by structural defects. Both in- 
frared and XRD measurements demonstrate that exter- 
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nal pressure can drive the MIT much closer to a genuine 
Mott transition, leading eventually to a homogeneous 
metallic phase, as expected insofar in (Vo.gsgCr . 011)203. 
However, to recover the lattice parameters and the elec- 
tronic structure of undoped V2O3 a much higher pres- 
sure is needed than expected from macroscopic transport 
measurements [H 0] ■ These findings shine new light on 
the properties of vanadium sesquioxide. In particular the 
phase separation occurring along some pathways across 
the PI-PM transition translates into intrinsically differ- 
ent properties of the corresponding PM phase. 



RESULTS 

In the PM state at 220 K, the optical conductivity 
cri(w) of (Vo.9 8 gCr . 011)203, reported in FigQJb) shows 
a pseudo-gap in the far infrared but absolutely no Drude- 
like contribution. One may also notice that the extrapo- 
lation to u—0 is in excellent agreement with the dc con- 
ductivity. The absence of a Drude term is in striking 
contrast with what is seen in V2O3 at the same tem- 
perature (Figfljc)), and implies that within the PM re- 
gion of the phase diagram the metallic properties can 
change dramatically. On the contrary, in the PI phase 
the spectrum of (Vo.gsgCro. 011)203 is akin to that of 
(Vo.g72Cro. 028)203: both spectra are characterized by a 
small insulating gap associated with an absorption band 
centered around 3000 cm -1 . The T-induced PI-PM tran- 
sition thus corresponds only to a weak red-shift of the 
spectral weight (FigQTb)). The optical data also show 
a strong thermal hysteresis (FigQTb)): u P on cooling the 
MIT occurs below ~ 260 K; upon heating thereafter the 
sample, the PM phase remains stable up to 275 K, and 
also at 300 K the PI phase is not fully recovered. 

The absence of a metallic contribution in the optical 
spectra of the x=0.011 sample calls for further investi- 
gations, especially since recent photoemission measure- 
ments on crystals from the same batch detected a clear 
metallic quasiparticle peak Tol- 13|. when performed with 
sufficient bulk sensitivity [l-jj. To better understand this 
apparently contradictory results, we performed SPEM 
with sub-micrometric spatial resolution [l5| , revealing an 
inhomogeneous electronic distribution in the PM phase. 
The maps of the photoelectron yield at Ep 
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are 

plotted in Fig. [2]together with representative photoemis- 
sion spectra. It should be pointed out that our experi- 
mental geometry [l5| and photon energy (27 eV) did not 
allow the collection of sufficient bulk sensitive photoelec- 
tron spectra Thus the spectra corresponding to 
the metallic phase do not show a pronounced quasiparti- 
cle peak. Nevertheless, the difference between insulating 
and metallic phases is unambiguous. Consequently, we 
interpret the photoemission microscopy images as repre- 
sentative of the topmost part of the metallic and insulat- 
ing domains, whose thickness is expected to be compa- 



rable to the IR probing depth. This penetration depth 
is sufficient to provide the excellent agreement between 
optical and the dc conductivity obtained with transport 
measurements. With these maps we resolve metallic and 
insulating domains which appear as T < XjyiiT and per- 
sist in the PM phase. A uniform insulating behavior is 
found instead in the PI phase (Figj2jc)), showing that 
the PM state only is strongly phase separated. Further- 
more, when the sample is cooled down again (Fig. [2fd)), 
the metallic and insulating domains are still found in the 
same position and with a similar shape, clearly showing 
that they form around specific points in the sample. It 
is natural to correlate these nucleation centers to struc- 
tural inhomogeneites in the material. Such a tendency 
towards phase separation is related to the correlation en- 
ergy which fosters a thermodynamic instability in the 
vicinity of the Mott transition [lj| [l9| .This is anologous 
to the liquid-vapour transition, and we can hence visu- 
alize the PI-PM transition as the formation of metallic 
clouds within an insulating sky. 

Although a T-induced structural phase separation was 
previously proposed on the basis of diffraction 0, 0| and 
EXAFS [5] experiments, it was never associated before 
with an inhomogeneous electronic state. Our data maps 
for the first time the MIT on the microscale, showing 
metallic and insulating domains which can be naturally 
linked to the a and /3 lattice structures, respectively. 

On the basis of these experimental findings in pho- 
toemission microscopy we can also understand the poor- 
metallic behavior in the PM phase, detected in the op- 
tical spectra. To further elaborate on this, we employ 
the Effective Medium Approximation (EM A) [2(| which is 
widely recognized to be meaningful in the presence of mi- 
croscopic inhomogeneities smaller than the wavelength of 
the electromagnetic radiation. EM A describes the (com- 
plex) d ielectric function se of an inhomogeneous system 
[2lL [22j in terms of the optical response of the two compo- 
nents (details in the supplementary material), weighted 
by the respective volume fraction / and depolarization 
factor q: 



f- 



q a 



+ (1-/)- 



£/9 - £E 



ep- 



q a 



Here, e a and ep are the complex optical dielectric func- 
tions of the metallic and the insulating phase, respec- 
tively. 

We reproduce the experimental conductivity in the PM 
phase of (Vo.gsgCro. 011)203 at any T by using as in- 
put parameters the a and (3 conductivities as measured 
in V 2 3 [iH and (Vo.972Cr .o28)20 3 , respectively. In 
particular previous measurements showed that the op- 
tical properties of V2O3 were strongly dependent on lat- 
tice modifications 23j. As the lattice parameters of a 
(Vo.gsgCro. 011)203 phase correspond to those of V2O3 at 
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400 K P, 13], we choose the corresponding 400 K opti- 
cal conductivity for the EMA reconstruction. Instead, 
the optical conductivity of (V .g72Cr . 028)203, represen- 
tative of the /3 phase, being nearly T-independent was 
choosen at 220 K. Extracting /= 0.45±0.05 from the 
map in FigJ^a) (see Supplementary Figure SI and dis- 
cussion in the supplementary material), and leaving q 
as the only free-fitting parameter a good agreement be- 
tween fit (empty squares in FigQJd)) and experimental 
data (solid line) is found for q= 0.30±0.05 at 220 K (q= 
0.40±0.05 at 300 K), similar to the value found in V0 2 
across its first order MIT j^. Since f>q the MIT tran- 
sition is percolative in nature, thus explaining the low 
dc values measured in samples coming from the same 
batch. Infrared measurements associated with SPEM 
data therefore clarify the poor-metallic behavior observed 
on a macroscopic spatial scale. 

To further support the paramount role played by 
phase separation for the low-energy excitations of 
(V0.989CT0. 011)203, we compare our infrared results with 
state-of-the art LDA+DMFT calculations H3] (Sup- 
plementary Figure S2 and discussion in the supplemen- 
tary material). The LDA+DMFT optical conductivity 
displays a metallic peak in the PM a phase (yellow line) 
and a gap followed by the t2 g — > <2g excitations in the 
insulating /? PI phase (gray line) (FigOTe)). To calcu- 
late the LDA+DMFT optical conductivity for the PM 
phase at x=0.011 we follow what our experimental anal- 
ysis suggests, namely to take a mixture of the pure el- 
and ^-phases. An important point is that, while EMA in 
the experimental analysis is (necessarily) performed by 
assuming the metallic phase to be represented by x=0 
and the insulating one by x=0.028, our LDA+DMFT 
approach allows us to use the "actual" a- and /3-phases 
at x=0.011 in the EMA formula. Using the same val- 
ues for the parameters / and q for theory and exper- 
iment, we find a LDA+DMFT optical conductivity at 
x=0.011 (open symbols in FigOJe)) which strongly re- 
sembles the experimental curve: in particular, its behav- 
ior is markedly insulating-like, as in the experiment. Note 
that tendencies towards phase separation close to a MIT 
have been discovered in DMFT before, on the basis of a 
negative compressibility [l8| and a negative curvature of 
the energy vs. volume curve These tendencies can 
be further enhanced by the coupling of electronic and 
lattice degrees of freedom. 

Let us now turn to the P-induced MIT: at the mini- 
mum P (1 kbar), cti(oj) shows an overall insulating be- 
havior resembling that at ambient pressure. As P is in- 
creased, 1 (cj) changes abruptly. Indeed, a sudden trans- 
fer of spectral weight from above ^3000 cm -1 to lower 
frequency builds up a metallic conductivity which well 
extrapolates to the dc values (solid diamonds in FigEJ) , 
measured on samples of the same batch. The difference 
from the poor-metallic behavior found in the T-induced 
PM phase at 220 K (and ambient-P) is striking. Never- 



theless, even at 15 kbar, the optical conductivity of V2O3 
is not yet recovered, challenging the proposed scaling re- 
lation between Cr-doping and pressure [x ~ 0.011 ks P~ 
4 kbar) obtained by dc resistivity measurements 0, H|. 
The metallicity induced by pressure can be associated 
with a double process: a weaker phase separation and a 
modification of the a-domain lattice parameters, which 
may become more similar to those of pure V2O3 at room 
temperature. 

In order to clarify these issues, high pressure XRD 
measurements have been performed in the same pressure 
range and the results for the (300) Bragg line (confirmed 
also by the similar behavior on other Bragg peaks) are 
shown in Fig|4ja). Two Bragg peaks appear within a 
range of ~ 1 kbar across the MIT, demonstrating that 
the PI phase coexists with the PM one only in this lim- 
ited pressure range. A further increase of P strongly 
reduces the phase separation (see FigHKb)) and above 10 
kbar the Bragg signal is made almost entirely by the PM 
component. 

These data indicate that a more homogenous elec- 
tronic state develops in the P-induced PM phase. Indeed, 
an EMA calculation reproduces the P-dependent optical 
conductivity (see inset of FigJ3]), through the same input 
curves as before. For the same value of g=0.4 obtained 
at ambient P and room T, a metallic volume fraction / 
=0.85 at 6 kbar and 0.95 at 15 kbar can be found. These 
values (solid symbols in FigSf b)) are in good agreement 
with values obtained from XRD measurements, but much 
larger than the metallic fraction / = 0.42 found at 220 K 
in the T-induced PM phase through SPEM. Eventually, a 
metallic volume fraction / ~1 is found at pressures where 
the c/a value in FigHJc) nearly attains that of pure V2O3 
(c/a ~2.83), thus suggesting the key role played by the 
lattice in the physical properties of V2O3 [231 l26j|. 



DISCUSSION 

We explored the phase diagram of (Vo.gsgCro. 011)203 
with a combination of different experimental methods 
and conjoined LDA+DMFT calculations, with particular 
attention to the PI-PM transition which can be obtained 
by changing doping, temperature or pressure. Our re- 
sults shine new light on the different metallicity of the 
PM phase obtained by lowering T with respect to the 
one obtained by applying an external P. We showed that 
a key role to explain this difference is played by phase 
separation: in particular, the coexistence of metallic and 
insulating domains was unambiguously demonstrated in 
the T-induced PM phase. After a thermal cycle these 
domains reappear in the same positions and with a sim- 
ilar shape: this memory effect may be related to nu- 
cleation centers, with electronic correlations further fos- 
tering phase separation. The PM phase reached under 
pressure is instead much more homogeneous and, by in- 
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creasing P, it soon covers a raising fraction of the volume. 
However, it is intrinsically different from pure V2O3 if one 
compares points which are nominally equivalent in the 
McWhan phase diagram, as demonstrated by the differ- 
ent lattice parameters. This suggests that, during the 
phase transition induced by P, the strain between the 
different lattice parameters between PM and PI seems to 
be accomodated in a more homogeneous way by reduc- 
ing their difference. In the T-induced transition, instead, 
this difference is unchanged and the strain is accomo- 
dated inhomogeneously inducing domains with different 
lattice and electronic properties. 

Our infrared results, obtained by changing doping, 
temperature and pressure, were captured correctly by 
employing the effective medium approximation. Thus, 
optical conductivity bridges between local probes like 
photoemission and XRD, which can directly detect the 
microscopic details of the phase separation, and trans- 
port measurements, which are instead sensitive to the oc- 
curance of a long scale metallic property. Thanks to this 
combination of different methods with intrinsically dif- 
ferent spatial scale sensitivity, we were able to clarify the 
properties of the phase diagram of vanadium sesquioxide 
from macroscopic to microscopic spatial scales, showing 
that the phase separation occuring along the different 
pathways across the PI-PM transition translates into dif- 
ferent properties of the PM phase. 

In conclusion, our investigation shows that the phase 
diagram of V2O3, celebrated for its essential simplicity, 
reveals on a microscopic scale a much richer phenomenol- 
ogy which has to be taken into account in the physical 
description of this model system. 



METHODS 
Materials 



extrapolation was performed by using a Drude-Lorentz 
fitting of the measured R by taking into account the 
dc values of the resistivity measured on samples of the 
same batch, while the high-frequency tail from Ref.[23[ 
was merged to the data, in order to perform reliable 
Kramers-Kronig transformations. A diamond anvil cell 
(DAC) equipped with type Ha diamonds was used to per- 
form pressure-dependent measurements. A hole of about 
300 /xm diameter was drilled in a 50 fim. thick Al gas- 
ket, allowing us to span over the desired pressure range 
(0-15 kbar) in a reproducible and cyclic way. A small 
piece was cut from the sample measured as a function 
of T and placed in the sample chamber together with a 
KBr pellet (used as hydrostatic medium) and the elec- 
tric field of the radiation was polarized within the ab 
plane. We took great care of creating a clean sample 
diamond interface. The pressure was measured in situ 
with the standard ruby fluorescence technique 27j. In 
order to reduce the experimental error of such data, a 
ruby micro-sphere was placed inside the sample cham- 
ber while a second one was placed on the external face 
of the diamond. Measurements were performed on both 
ruby balls several times for each pressure, taking there- 
fore in account any slight change in the temperature or in 
the spectrometer calibration. The incident and reflected 
light were focused and collected with a cassegrain-based 
optical microscope equipped with a MCT detector and 
a bolometer and coupled to a Michelson interferometer. 
This allowed us to explore the 350-12000 cm~' spectral 
range exploiting the high brilliance of Synchrotron Radi- 
ation at SISSI beamline at ELETTRA storage ring [28|. 
The measurement procedure was the same as described 
in Refs. [29|, [3(|. The conductivity, cr(w), has been cal- 
culated through Kramers-Kronig transformations, taking 
into account the diamond-sample interface 3ll433| . Data 
are not shown in the 1600-2700 cm -1 due to the strong 
phonon absorption of diamonds. 



All the measurements on (V0.989CT0. 011)203, employ- 
ing different experimental techniques as detailed in this 
paper, were performed on single crystals from the same 
batch grown at Purdue University. The specimens were 
thoroughly characterized prior to the experiments by 
means of transport measurements (Patricia Metcalf, pri- 
vate communication), and precisely oriented with X-ray 
diffraction. 



Infrared Measurements 

The temperature-dependent reflectivity (R) measure- 
ments have been performed at near-normal incidence 
with the electrical field polarized within the ab plane 
and by using as a reference a metallic film (gold or alu- 
minium depending on the spectral range) evaporated in 
situ over the sample surface. A standard low-frequency 



Photoemission Microscopy 

The photoemission microscopy results were obtained 
on the Spectromicroscopy beamline [l5] at the Elettra. 
Photons at 27 eV were focused through a Schwarzschild 
objective, to obtain a submicron size spot. The 
measurements were performed in ultra-high vacuum 
(P< 2xI0~ 10 mbar) on in situ prepared surfaces of 
(Vo.989Cr . 011)203. Photoelectrons were collected at an 
emission angle of 65 degrees from the normal, which 
makes the photoemission spectra and images surface sen- 
sitive. A standard photoemission microscopy procedure 
was employed to remove topographic features j3~i| from 
the images presented in FigfJl and to obtain an unam- 
biguous contrast between metallic and insulating regions 
[l6| : the spectral intensity at the Fermi level ~Ep - high in 
the PM and absent in the PI phase ll|, [l2j - was divided 
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by the intensity recorded at the V3d-band. Since the lat- 
ter remains essentially constant between the PI and PM 
phases, two-dimentional maps of this ratio give a genuine 
representation of the lateral variation of metallic and in- 
sulating domains. 



X-ray diffraction 

The X-ray diffraction measurements were performed 
using 17 keV photons dispersed by the CRISTAL beam- 
line at the SOLEIL synchrotron radiation source [35| . 
Pressure was applied to a (Vo.gsgCro. 011)203 powder by 
means of a membrane driven diamond anvil cell, and its 
value in the cell was carefully measured during the mea- 
surement sequence using the standard ruby fluorescence 
method. The transmission of the diamond cell at the 
photon energy employed during these experiments was 
sufficient to allow a high resolution detection of the main 
Bragg reflections like the (300) and (214) presented here. 
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Figure 1: Temperature dependent infrared measurements, (a) Schematic phase diagram as defined in Ref. 0]; markers 
indicate the measurements displayed in (b)-(e) in the same color coding, b) <ji(lu) of (Vi-zCr^Os for x=0.011 and 0.028 
at selected T. cti(cj) for x=0.011 displays different behaviors in the PI and PM regions. The two curves measured at 300 
K highlight the presence of large hysteresis: the curve at 300 K reached from the PI side (upon cooling) shows the expected 
insulating behavior that is not fully recovered when 300 K is reached from the PM phase (upon heating). In the PI phase only a 
small difference is found in ai(uu) of x=0.011 and 0.028. Solid symbols in panels b to d represent dc values as measured in Ref.Q 
on crystals of the same batch, (c) g\(uj) in the PM phase highlighting the difference between V2O3 and (Vo.gggCro. 011)203 
showing a clear Drude peak and a poor-metal behavior, respectively, (d) Experimental and (e) theoretical reconstruction 
of <7i(w) in the T-induced PM phase, succesfully obtained by modelling, with the same q and / parameters (see text), the 
coexistence of metallic and insulating domains through the effective medium approximation (EMA). 
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Figure 2: Temperature dependent photoemission microscopy measurements. Scanning photoemission microscopy 
images and spectra on (Vo.989Cro. 011)203, collected at 27 eV photon energy and at different temperatures on a 50 /j,m by 
50 /jm sample area. The images are obtained with standard procedures to remove background and topography effects [fa ]. 
The pictorial contrast between metallic and insulating zones is obtained from the photoemission intensity at the Fermi level 
[f7| . normalized by the intensity on the V3d band (binding energy -f.2 eV). The lateral variation of this normalized intensity 
is represented in the color scale (the intensity scale is the same for all four maps) from its minimum value, in violet, to its 
maximum, in red, and is therefore a direct visualization in real space of the metallicity of the system, fnhomogeneous properties 
are found within the PM phase at T=220 K (a) and 260 K (b), where metallic (in red) and insulating (in blue) domains coexist. 
To recover the PI phase, the sample was heated to 320 K, to make sure that the hysteresis effects still present at 300 K (see 
Fig. 1) are overcome: at this temperature (c) a homogenous insulating state is obtained. After a whole thermal cycle (d) the 
metallic regions can be found in the same position and shape as in (a). The photoemission spectra (outer panels) from selected 
representative areas (A,B,C) corroborate this interpretation. 
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Figure 3: Pressure dependent infrared measurements. Optical conductivity of (Vo.gsgCro. 011)203 under pressure (data 
between 1600 and 2700 cm -1 , due to the strong diamond phonon absorption, are substituted by markers as a guide for the 
eye). Under increasing P, the optical conductivity is depleted at energies above 3000 cm -1 and increases below this energy. 
The curves at 1 kbar, 6 kbar and 15 kbar are compared to the PM state of V2O3 at 300 K and to the poor metallic phase 
of (Vo.9S9Cro. 011)203 at 220 K. In the latter phase no metallic-like increase is built up for oj — > 0, whereas as a function of 
pressure some coherent transport takes place, as supported also from the resistivity data (solid diamonds in Fig(3| that allows 
to extrapolate smoothly the curves to uj — s> 0. In the inset is shown the experimental reconstruction of eri(w) in the P-induced 
PM phase obtained by EMA. 
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Figure 4: Pressure dependent Bragg X-ray diffraction measurements, (a) Bragg X-ray diffraction (300) line of 
(Vo.9S9Cro. 011)203 under pressure at room temperature. The peak at the left corresponds to the lattice parameters of the PI 
phase, while the one at the right corresponds to PM. The spectra from violet to red were taken at pressures from to 20 kbar, 
corresponding to the data points in Fig(4jb) and (c). (b) PM intensity over the total one (PM+PI) as a function of P for the 
(300) line (circles) and for the (214) line (triangles). Solid symbols represent the same quantity calculated from the / EMA 
parameter (see above), (c) Ratio of the lattice parameter c/a vs. P extracted from Figfija). 
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Supplementary material 

A Microscopic View on the Mott tran- 
sition in Chromium-doped V2O3 



S. Lupi, L. Baldassarre, B. Mansart, A. Perucchi, A. Bari- 
nov, P. Dudin, E. Papalazarou, F. Rodolakis, J. -P. Rueff, 
J.-P. Itie, S. Ravy, D. Nicoletti, P. Postorino, P. Hansmann, 
N. Parragh, A. Toschi, T. Saha-Dasgupta, O. K. Andersen, 
G. Sangiovanni, K. Held and M. Marsi 

Effective medium approximation 



Within the Bruggeman EMA theory (Bruggeman, D.A.G., 
Ann. Phys. Lpz., 24, 636 (1935)), which has been success- 
fully used to describe the optical properties of phase-separated 
correlated systems (Ref.22, main paper), the complex optical 
dielectric function (ee) is described by the following equation: 



' EE 



+ 



+ (1-/) 



£/3 — £E 

£0 + ^f-£E 







where l a (/) and ep (1 — /) are the complex optical dielectric 
functions (the filling factors) of the metallic and insulating 
phase respectively and q is the depolarization factor that de- 
pends on the shape of the metallic domains dispersed in the 
insulating matrix. 

We reproduce the experimental conductivity in the PM 
phase of (Vo.989Cro. 011)203 at any T by using as input pa- 
rameters the a and /3 conductivities as measured in V2O3 
(Ref. 23, main paper) and (Vo.972Cro. 028)203 (Fig. lb in the 
paper), respectively. 

The filling factor was determined from the photoemission 
images of the spectral weight at Ef (Fig. 2 in the paper). In 



particular the bimodal histogram of the intensity distribution 
extracted from the PM map (Supplementary Fig.Sl) shows 
unambiguously that at 220 K two main phases can be iden- 
tified: one insulating, corresponding to the violet-blue colors 
in the photoemission map, and one metallic, corresponding 
to orange-red colors. The broad minimum between these two 
components corresponds to the green color in the microscopic 
maps, which we interpret as regions where metallic and insu- 
lating regions are too small to be resolved with our instrumen- 
tal resolution. These considerations lead us to conclude that 
the assumption of a two-component description for the PM 
phase is perfectly legitimate. Furthermore, the PM histogram 
can be fitted by a sum of two Gaussian curves (dashed line), 
whose relative intensities provide the insulating and metallic 
filling factors. Specifically, the fit provides for the metallic 
part f=0.45±0.05. Remarkably, at 320 K only (not shown) 
one Gaussian (corresponding to the insulating phase) is suf- 
ficient to fit the intensity distribution histogram, as expected 
from the corresponding photoemission microscopy map (Fig. 
2c in the paper), which shows an almost uniform blue- violet 
color. 

As the EMA analysis is strongly constraint by the 
pseudogap behavior of the optical conductivity below 1000 
cm -1 (see Fig. Id in the paper) as well as by the finite 
(independently measured) dc value (cr(dc) ~ 250 SI -1 cm -1 ), 
by changing / in the EMA analysis accordly to the previous 
estimate, the reconstruction is still possible with a similar 
quality for different values of q (as far / > q in order to have a 
finite dc conductivity). From f=0.45±0.05 one can estimate 
q=0.30±0.05. 



LDA+DMFT Optical conductivity 



The formula for the optical conductivity reads 



J 



5Ra a/3 (n) = 



2ne 2 h 
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,f(L j')-f(,w'+Q) 
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where V is the volume of the unit cell, a and ft are the 
Cartesian coordinates x, y, z, f is the Fermi distribution 
and finally v (k) denote the Fermi velocities obtained from 
the LDA Hamiltonian e LDA (k) . The spectral function A reads 



A_(k, w) = - -Im [(w + n)l - | LDA (k) - £(w) 



(2) 



Here, the double underlines expresses the matrix character 
in spin and orbital space. Hence, the correlation effects 
enter Eq. H) only via the DMFT self energy E(oj) in A(k)- 
Let us also note that due to the trace in Eq. fl]) orbital off 
diagonal contributions are fully taken into account. Finally 
we remark that in expression ^ vertex corrections are 
neglected. For the single band case, this becomes exact in 
the limit of infinite dimensions, just like DMFT does. In a 
multi-band system with similar orbitals, such as V2O3 it is 



a good approximation (Tomczak, J. M., and Biermann, S., 
Physical Review B 80, 085127 (2009)). 

Calculation of S and A 



Our LDA+DMFT calculations were carried out on the 
crystal structure of Vo.9S9Cro. 011)203 in the a and j3 phases 
( Robinson, W. R., Acta Crystallographica Section B, 31, 
1153 (1975) ). The LDA calculations were performed in the 
localized basis of muffin tin orbital (MTO) based LMTO and 
NMTO ( Andersen, O. K. and Jepsen, O., Phys. Rev. Lett., 
53, 2571 (1984); Andersen, O. K. Lecture Notes in Physics 
Springer, New York, 2000)). The full LDA Hamiltonian was 
downfolded on to the V t2g only low-energy Hamiltonian, 
defined in the effective Wannier function representation, 
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using the NMTO-downfolding technique. This yields a 12 
x 12 Hamiltonian for each k-point of the Brillouin Zone as 
there are four vanadium atoms in the primitive unit cell. For 
a detailed discussion of the downfolding procedure and the 
setup of the LDA+DMFT calculation (Saha-Dasgupta, T., et 
al, arXtv.org, 0, 0907.2841 (2009); Ref.26, main paper). For 
further details about LDA+DMFT see Held, K., Advances 
in Physics, 56, 829 (2007). To obtain the self energy 
matrix S(cj) wc first perform an analytical continuation 
of the local Green function G loc '(o;) using the maximum 
entropy code by Jarrell and Gubernatis (Jarrell, M. and 
Gubernatis, J. E., Physics Reports, 269(3), 133 (1996) ) 
and some improvements described by Gunnarson, Haverkort, 
and Sangiovanni (Gunnarsson, O., Haverkort, M.W., and 
Sangiovanni, G.,Phys. Rev. B 81, 155107 (2010)). In a 
second step we solve a multidimensional rootfinding problem 

^ (M) = ^// fc ( M + ,)l-^)-£M (3) 



In Supplementary Fig.S2 we report the orbital-resolved 
LDA+DMFT self energies. In the plots we set the Fermi 
energy to £f = and plot the self energy in green and 
the ai g self energy in blue. We summarize the quantities for 
the a- and the /3-phase on the left hand and right hand side 
of the panels respectively. Overall our results agree with the 
results of the previous LDA+DMFT analysis (Ref. 26, main 
paper), although we performed the calculations at slightly 
lower U = 4.0 eV values (see Ref.26, U = 4.2 eV) in order 
to get the best agreement with experiment. Indeed we found 
an extreme sensitivity of the optical gap depending on the 
value of the Hubbard U. For a detailed analysis of the U de- 
pendence of the LDA+DMFT calculations see Toschi, A., et 
al, Journal of Physics: Conference Series, 200, 012208 (4pp) 
(2010). 
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Supplementary Fig.S 5: Intensity distribution of PM 
photemission map. Histogram of the intensity distribution 
of the PM photoemission map (ratio between the photoemis- 
sion intensities at E_f and on the V3d band) (open symbols) 
as extracted from the map at 220 K in Fig. 2a in the paper. A 
bimodal distribution can be observed. This distribution was 
fitted with two gaussian components (dashed line) providing 
f=0.45±0.05. The color scale on the horizontal axis is the 
same used in the photoemission maps in Fig. 2 of the paper. 
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Supplementary Fig.S 6: LDA+DMFT self energy. 

LDA+DMFT self energy for (V .989Cr .oii)203 in the a 
phase (left hand side) and the /? phase (right hand side) for 
U — 4.0 eV and T — 390A". In the top and bottom panels we 
show the LDA+DMFT self energy for the and ai g states 
respectively. 



